Role of spin-orbit coupling in the far infrared absorption of lateral semiconductor dots 
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We investigate the relevance of the spin-orbit coupling to 
the far-infrared absorption of two-dimensional semiconductor 
dots. Varying the strength of the Dresselhaus term, a mech- 
anism feasible in experiment by changing the dot width, dis- 
tinctive splittings of the Kohn peak as well as additional low 
energy modes are predicted in a non-interacting model. Each 
mode has a spatial distribution of charge, perpendicular and 
in-plane spin densities that correlate in a peculiar way with 
the frequency and polarization of the external field. We study 
the robustness of these features against electron-electron in- 
teractions as well as the appearance of interaction-induced 
additional characteristics. 

PACS 73. 21. La, 73.21.-b 



The role of the electron spin-orbit (SO) coupling in the 
general properties of lateral semiconductor dots has at- 
tracted much interest in recent years. To a large extent 
this is motivated by the technological impact attributed 
to the so-called spintronic devices, whose operation prin- 
ciple exploits precisely the SO mechanism to control the 
spin transport jlj. Actually, spintronics is an incipient 
field whose fundamental principles need a deeper under- 
standing In this general context, experimental and 
theoretical investigations have very recently shown that 
the SO coupling affects the quantum transport and, more 
specifically, the conductance fluctuations of chaotic quan- 
tum dots in a parallel magnetic field ^-^|. 

The different SO contributions in semiconductor dots 
have been reviewed by Voskoboynikov et al. In gen- 
eral, SO coupling has two contributions: the Rashbaterm 
due to the nanostructure confinement potential; and the 
Dresselhaus term originating from the crystalline inver- 
sion asymmetry. In lateral dots the strength of the latter 
depends on the vertical width, which can be experimen- 
tally varied with vertical electric fields. For an embed- 
ding GaAs medium both Rashba and Dresselhaus param- 
eters are reasonably determined 0. 

In this Letter we study the influence of SO couplings 
on the dipole absorption of GaAs dots with a parabolic 
confinement. While the Rashba term provides a negligi- 
ble effect we show that the Dresselhaus contribution leads 
to characteristic splittings of the Kohn modes and to the 
appearance of quasi-spin flip modes in the far-infrared 
absorption (FIR). For each mode we find a correlation 
between the spatial patterns of oscillating electron den- 
sity, horizontal (in-plane) spin density, and vertical spin 
density. As in most of the literature treating the SO 
coupling j|-U, we shall first consider a non- interacting 



model, using both analytical and numerical treatments. 
Nevertheless, the role of the interaction will be also con- 
sidered within self-consistent mean-field for Hartree and 
exchange-correlation contributions, using a spinorial for- 
malism of local-spin-density theory (LSDA) || . 

Within the effective Hamiltonian formalism and for a 
quasi two-dimensional system the Rashba and Dressel- 
haus [for the standard (001) plane of GaAs] contributions 
read, respectively, 
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where the A's control the coupling strengths, the er's re- 
fer to the Pauli matrices and V(r) stands for the to- 
tal electrostatic potential. In Eq. P gives the ki- 
netic momentum. Although we shall not consider in this 
work an external magnetic field (or vector potential A) 
it can be easily included using the canonical momentum 
P = -ihV+ f A in Eq. (0). 

For GaAs one has A# ; = 5.33 A 2 and a material-specific 
constant 7 = 27.5 eVA 3 The latter determines the 
Dresselhaus coefficient in terms of the dot vertical width 
zq as [|| Xd ~ j(tt/zo) 2 ■ As a representative case, we 
consider an external parabolic confinement with energy 
TujJq = 4.2 meV and a number of electrons ranging from 
6 to 12. A strict two-dimensional motion is assumed 
and the dot width is simply taken as a parameter for 
the Dresselhaus coupling. The Rashba term turns out to 
be negligible in all cases we have studied. However, the 
Dresselhaus term increases its importance for decreasing 
doth widths, leading to relevant effects as we shall show 
below already for widths zq w 50 A. The influence of 
the Rashba term has already been studied in the liter- 
ature for InSb dots 0, where band non-parabolicity is 
important, and in phenomenologically fitted models pM . 

We have implemented the numerical solution on a spa- 
tial grid of the full Hamiltonian TCfuii = 7~to + T~Cr + Ti-D, 
with 
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However, before presenting numerical results it is worth 
to obtain an approximate analytical solution of the rel- 
evant contributions 7io + T~Cd = "H with an approach 
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similar to that of Aleiner and Fal'ko H. A diagonaliza- 
tion to second order in Ad is obtained by means of the 
transformation 7i = U + HU, where 



U = cxp 
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We find 
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Identifying the orbital angular momentum £ z — xP y — 
yP x in (0), we realize that the Dresselhaus term intro- 
duces a correction depending on l z a z quite similar to 
that of the Rashba term for a circular system, but in the 
transformed reference frame. In fact, in the laboratory 
frame the Ti eigenstates to first order in are spinors 
that slightly deviate from the a z eigenspinors. We shall 
call them quasi-up Xqu and quasi-down Xqd states 



Xqu = <Pnt{?) 
Xqd = ^nf( r ) 
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where (p n e(r) are the eigenfunctions of the circular har- 
monic oscillator H.q. Since the dipole operator is invari- 
ant by the transformation (|^), the analytical solution pre- 
dicts dipole transitions with A£ — ±1 for quasi-up and 
quasi-down subsystems (see Fig. 1) at the energies 
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An inmediate consequence is the splitting of the original 
Kohn mode at hujQ in two peaks separated by an energy 
Ahio = 2X 2 D m/h 2 . 

Left panels of Fig. 2 show the numerical result for the 
FIR absorption of TLf u ii obtained within a real-time sim- 
ulation after an initial rigid displacement of the system 
as in Ref. jnj. The results correspond to N = 6 and 10 
electrons with a Dresselhaus parameter for a dot width 
of zq = 50 A. The splitting around Tiojq nicely agrees 
with the analytical prediction (arrows), thus proving the 
validity of the model. For the N — 10 case we note the 
existence at low energy of a peak with much less strength 
than the dominant ones (notice the enlarged scale). This 
too can be explained within the analytical model realiz- 
ing that N = 10 has a partial occupation of the third 
oscillator major shell and thus the possibility of quasi- 
spin flip transitions \qu Xqd at an energy 4:X 2 D m/h 2 , 
indicated with dotted lines in Fig. 1 and with an arrow 



in Fig. 2. Note that in spite of having A£ z = these 
transitions can manifest in the dipole spectrum through 
the higher order corrections in Eq. (||). 

Due to the spin-orbit coupling the charge dipole oscil- 
lation also induces a spin oscillation. This is proved in 
the right panels of Fig. 2 where the Fourier transform of 
(S x + S y ) is plotted. The same frequencies of the dipole 
absorption manifest in the spin channels, although with 
different relative strengths. It is also interesting to ana- 
lyze the spatial patterns of induced density Sp(r,t), par- 
allel spin 5Sii(r, t) and vertical spin SS z (r,t) for a given 
oscillation mode and a fixed time (Fig. 3), where we de- 
fine the spin density as S(r) = £\ =1 (xi | <?<5(r - n) | Xi)- 
This local-signal analysis has been performed using the 
method of Ref. |T|] . Quite remarkably, the oscillation of 
6S z (r,t) is always in a perpendicular direction to that of 
the charge density, which is explained in the analytical 
model as a coherent quasi-up and quasi-down excitation 
due to the SO term. On the contrary, the parallel spin 
patterns depend on whether the mode has quasi-spin flip 
character or not. For the split Kohn modes (peaks a and 
b) <JSy(r, f) is localized in the two regions perpendicular 
to the charge oscillation axis. However, the quasi-spin 
flip peak (c) shows no angular dependence of the pat- 
tern. It can be shown that this behaviour is in complete 
agreement with the prediction of the analytical model 
using the spinors of Eq. (||) . 

We consider next how the features discussed above are 
modified by the addition of the Coulomb interaction 
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where k is the semiconductor dielectric constant (k = 
12.4 for GaAs). We shall estimate interaction effects us- 
ing a selfconsistent mean-field for the Hartree electro- 
static potential 
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and the LSDA for electronic exchange and correlation. 
Within the present spinor formalism the dots have a vec- 
torial spin magnetization S(r) and, therefore, a general 
exchange-correlation functional of the type Exc[Pi S] is 
necessary. For this we have resorted to the non-collinear- 
spin formulation of LSDA theory, described for instance 
in Ref. ||, an d employed the Monte-Carlo results of Ref. 
]T3] | for the bulk electron gas. 

The modifications stemming from the Coulomb in- 
teraction are twofold; first, the electron-hole transitions 
move in energy due to the change in the effective mean 
field and; second, new collective peaks may appear as 
a result of coherent electron-hole excitations. Figure 4 
shows the interacting results for N = 6 with zq = 50 
and 35 A. Note that the splitting of the Kohn-like modes 
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disappears for the larger width and, at the same time, 
a collective peak begins to appear at an energy below 
the electron-hole states. These results clearly show the 
formation of the collective plasmon and magnon states 
lying above and below the electron- hole transitions, re- 
spectively. We stress that in the absence of SO coupling 
the magnon state does not contribute to the FIR ab- 
sorption. The lower left panel of Fig. 4 indicates that 
a great enhancement of the magnon contribution is ob- 
tained by reducing the dot width, which is accompanied 
by an important Landau damping of the plasmon state 
into a bundle of peaks. The Coulomb interaction thus 
suppresses the splitting of the Kohn-like modes although 
it introduces the mechanism of Landau damping that ul- 
timately, and rather abruptly with zq, dominates for de- 
creasing widths. 

The right panels of Fig. 4 show the horizontal spin 
spectra induced by the dipole shift. These results permit 
to ascertain the spin character of the different modes and, 
by comparison with the left panels, quantify their rele- 
vance to the FIR absorption. The relative spin or density 
character of each state is also seen from the amplitudes 
A of density and spin oscillation. For instance, when 
zq = 50 A we find A[8S z {v,t)] w 60A[Sp(r,t)] for the 
magnon and A[SS z (r,t)] w 0AA[Sp(r,t)) for the plas- 
mon, which shows the dominant spin and density type 
for magnon and plasmon states, respectively. 

The characterization of the interacting peaks is facili- 
tated by the comparison of the spatial patterns with the 
corresponding ones in the non-interacting model. Figure 
5 displays the density and spin patterns for two of the 
peaks of the zq = 35 A spectrum (lower left panel of Fig. 
4). As in the non-interacting model, the oscillation of 
SS z (r, t) is orthogonal to that of Sp(r, t) and the horizon- 
tal spin patterns show a priviledged spin direction, i.e., 
they have a net induced spin. Note in particular that 
peak b of Fig. 4 has no angular dependence in the hori- 
zontal spin pattern while for peak a there are two regions 
with enhanced 8S\\ (r, t). We have checked that these fea- 
tures are not changed when the polarization of the mode, 
given by the direction of the density oscillation, is varied. 
We stress that the mechanism discussed here would per- 
mit a control of the horizontal spin oscillation by tuning 
the frequency of the applied electric field. This effect of 
spin-to-charge conversion is a major requisite for hybrid 
spintronic devices B. 

Finally, Fig. 6 presents the spectra with Coulomb inter- 
action obtained for N = 10, an open-shell system in the 
analytical model. At zq = 50 A the open-shell character 
of this dot is clearly maintained, with low-energy elec- 
tron hole transitions and an associated collective mode 
at lo Ri 0.7 meV similar to the quasi-spin flip excitations 
of Fig. 1. For narrower dots, however, the energy of quasi- 
spin flip transitions increases and as a result this mode 
merges with the magnon and other electron-hole states. 
Contrary to the quasi-spin flip state, both magnon and 



plasmon mean energies are not dependent on the Dres- 
selhaus coupling. 

To summarize, the role of the SO coupling in the FIR 
absorption of GaAs dots has been studied. In an ana- 
lytical model without interaction the relevant Dressel- 
haus contribution induces characteristic splittings and 
quasi-spin flip modes. The Coulomb interaction favors 
the formation of collective plasmon and magnon states 
and it suppresses the splittings of the analytical model. 
Nevertheless, below a certain dot width Landau damping 
dominates the spectrum. Each mode has characteristic 
patterns of induced charge, horizontal spin and vertical 
spin, as well as a ratio of oscillating density and spin 
amplitudes. Therefore, by tuning the frequency of the 
applied electric field one can effectively control the dy- 
namical spin properties and the relative intensities of the 
charge and spin oscillation. 
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quasi-up quasi-down 
FIG. 1. Energy level scheme within the analytical model 
for Ho + Hd- The Fermi level for N = 6 is indicated by e F 
and the vertical arrows correspond to the allowed transitions 
having Al = ±1 within each quasi-spin subset. The dotted 
lines show quasi-spin flip Al = transitions within the N = 2 
major shell, discussed in the text for N = 10. 
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/i<2)(meV) hm(me\l) 
FIG. 2. FIR absorption (left panels) obtained from the 
Fourier transform of the real time dipole oscillation with TLfuii 
for a dot width of zo — 50 A. Right panels display the cor- 
responding frequency transform of the horizontal spin signal. 
The arrows indicate the energies of the split Kohn modes (at 
~ 4 meV) and of the quasi-spin flip modes (at w 0.5 meV) 
using the analytical model for this zq. 

FIG. 3. Spatial patterns of oscillating density 8p(r,t), 
horizontal spin <5Si|(r, t) and vertical spin SS z (r, t) at a fixed 
time for the a and c dipole modes of Fig. 1. The gray color 
scale of left and right panels indicates the signal magnitude 
while the corresponding sign is superimposed in white. 
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FIG. 4. Same as Fig. 2 including Coulomb interactions. 
The dashed line in the left panels displays the electron-hole 
transitions in the selfconsistent mean-field. 



FIG. 5. Same as Fig. 3 including Coulomb interaction. 
Upper and lower row correspond to the peaks labeled a and 
b in Fig. 4, respectively. 
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FIG. 6. Same as Fig. 4 for N = 10. 
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